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Abstract. We explicitly determine the large deviation function of the energy 
flow of a Brownian particle coupled to two heat baths at different temperatures. 
This toy model, initially introduced by Derrida and Brunet 1 , allows not only to 
sort out the influence of initial conditions on large deviation functions but also to 
pinpoint various restrictions bearing upon the range of validity of the Fluctuation 
Relation. 
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The recent upsurge of interest in nonequilibrium statistical mechanics lies in the 
discovery of simple yet generic results, embodied by the Fluctuation Relation initially 
brought forth by Evans, Cohen and Morriss [2], and later formalized into a theorem by 
Gallavotti and Cohen [3] . This Fluctuation Relation states that in the nonequilibrium 
steady-state of a (chaotic, Anosov) dynamical system, the temporal large deviation 
function of an appropriately defined entropy current verifies a particular symmetry 
property under time reversal. In practice the latter entropy current consists of a heat 
or particle flux and could be accessed experimentally in some situations @j. Before 
going into the specifics of that relation, we mention that a Markov dynamics analog 
of the Gallavotti-Cohen theorem was then found by Kurchan 5. and Lebowitz and 
Spohn 6 . More recently, other contributions on the subject have been put forward 
by Van Zon and Cohen and by Bonetto et al. J^j, where an extension of the 
Fluctuation Relation to quantities other than the entropy flux is proposed. 
The purpose of the present report is to present a simple model for which 
nonequilibrium fluctuations can be investigated precisely at finite times. We will 
show that a seemingly innocuous initial condition dependence restricts the domain 
of validity of the Fluctuation Relation, even though the latter applies to asymptotic 
times. 

In order to illustrate our point, we exploit a toy model introduced by Derrida and 
Brunct in a recent pedagogical account £Q. After having described this model, we 
will recall the standard statement of the Fluctuation Relation, and then we shall show 
that some initial condition effects, that intuition suggests to discard, can lead to an 
actual failure of the Fluctuation Relation. 

We take a particle in contact with two heat baths imposing distinct temperatures T\ 
and T2, whose velocity v(t) evolves according to the Langevin equation 

row = -(7i+7a)w + £i(i)+&(t) , (1) 
with =0 and 

= 27i2iM(* - = 2!><M(* ~ . ( 2 ) 
where the Boltzmann's constant has been taken equal to one. Such an equation can 
also describe the dynamics of apparently different models, for example, the dynamics 
of a hard rod connected at both ends to two heat baths (as a slight variation of 
the model introduced in [Sj , and is sketched in Fig . We also note that another 
model, consisting of two overdamped brownian particles interacting with a harmonic 
potential can be exactly mapped on Eq. Q). Finally, we quote that other 

analogies in term of electric circuits, in the spirit of J2| can a ls° be described by a 
Stochastic differential equation of the same type. Strictly speaking, the dynamics of 
the particle is an equilibrium one, since in Eq. the sum of the Gaussian noises is 
another Gaussian noise with viscosity 7 = 71 + 72 and noise strength D = Di + D 2 - 
The velocity probability distribution function (pdf ) of the particle is hence a Gaussian 
with zero- mean and variance T = D/j. Nevertheless, for our purpose, we distinguish 
the thermostats and we focus on the heat flowing from one thermostat (say for example 
thermostat 1) to the particle over a time duration t. This quantity is exactly the total 
work performed by the thermostat, and it reads: 

Qi= f dr«(r)(ei(r)-7i«(T)) , (3) 
Jo 

where i = 1,2. The pdf of this integrated injected power P(Qi,t) is the central object 
of our note. The Fluctuation Relation, applied without care, states without restriction 
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Figure 1. (color online) Sketch of a model, inspired by 9 described by the 
Langcvin equation The two particles are constrained to move together along 
the horizontal axis, and are in two heat baths at different temperatures. 



that for the fluctuating time averaged flux q — Qi/t, the large deviation function (ldf) 
ir(q) defined as 

n(q) = lim -\nP(qt,t) (4) 

f^oo t 

verifies the symmetry property 

- n(-q) = ■ (5) 

This is the Fluctuation Relation (FR). In Eq. ©, the constant e = ^- — tjt plays the 
role of an external field driving the heat flux. In the following, we will also consider a 
restriction of the relation JSJ to a finite interval (i.e. for \q\ < q* , where q* is a positive 
constant), in analogy with deterministic systems |13ll%| and stochastic systems [7]. We 
will refer to this last relation as "Extended Fluctuation Relation" (EFR), from the 
terminology of 0. Usually the Fluctuation Relation is supposed to be verified by 
the entropy flux S = Q\/T\ + Q2/T2, which differs from eQi only by a total time 
difference, and hence, following |HJ, Qi should, at least, verify an Extended Fluctuation 
Relation. 

We now wish to explicitly determine P(Qi,t) at all time and investigate the 
properties of the large deviation function. The techniques used are similar to those 
employed by Farago [T3|. We first introduce the joint probability p(v,Q\,t) that the 
particle has a velocity v and a given value of Qi at a time t. The pdf P(Qi,i) is 
obtained by integrating p over the velocity v. We also define the Laplace Transform 

p(v,X,t) = J dQ 1 e- x ^p(v,Q 1 ,t) , (6) 

which verifies: (i) p(v,0,t) = f(v,t) with f(v,t) the velocity pdf of the brownian 
particle, and (ii) P(X,t) = J dvp(v, A, t), with P(X,t) the Laplace Transform of 
P(Qi,t). The pdf p(v, A, t) verifies a Fokker-Planck Equation (FPE): 

d 

— p[v, A, t) = L x p(v, A, t) , (7) 
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In order to get the solution of the above FPE we first note that the eigenvalue equation 

L X f n (v,\)=Hn(\)fn(v,\) (9) 



is solved by: 



f..W = |(l-(l + 2n>,), (10) 



where: 



/ 4A 

r) = Jl + —(j 2 D 1 -j 1 D 2 -XD 1 D 2 ), (12) 
2D 

r ( A > = 7 (l + ,) + 2AZV T *^= T ^> (13) 

and where H n (x) denotes the Hermite polynomial of order n. We remark that the 
largest eigenvalue no(X) presents the symmetry 

Mo(A) = Mo(e-A) , (14) 

with e = 1/Tg — 1/Ti- Should 7r(g) be the Legendre Transform of ^o(A), this last 
relation would exactly yield |J5J . The solution of (JJJ for a given initial condition is 
hence easily obtained as: 

00 

p(v, \, t\v ) = J2 e^ Wt C n (v, X\v Q )f n (v, A) , (15) 

n=0 

where C n (v, X\vq) is the projection of the n-th eigcnfunction onto the initial state, 
which has been chosen to be a Dirac function centered in v = vq. In order to simplify 
the notations we will set all the scales to unity (i.e. 7 = D = Tq = 1), and introduce 

71 I + A7 . 72 I-A7 

7i = — = o ' 72 = — = 7, . ( 16 ) 

72 72 

- Di 1 + AD ~ D 2 1 - AD 

where both A7 and AD take values between —1 and 1. Thus, when T\ = T 2 one has 
that A7 = AD. Finally, integrating p over the velocity, one obtains the expression of 
P(X,t) for a given initial velocity vq: 

s ' ( w x /w » ™ ,sinh(?rf)V* / A(l + A 7 + (1 + AD)A) \ 

P ( A ' ^) = ^ + + + ) 6XP (f ( l + A(l + AD) +7? coth(^) J 

\ A detailed derivation of this operator for a slightly different system is given in: P. Visco, A. Puglisi, 
A. Barrat, E. Trizac and F. van Wijland, submitted to J. Stat. Phys, cond-mat/0509487 
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The long time behavior of P(X, t) is clearly dominated by P(X, t) ~ e M °( A )*, with /Zo(A) 
the largest eigenvalue of L\. This result was already found in pQ. The expression of 
Ho presents two cuts in the real axis for A > A + and A < A_, with: 



AD-Aj 1-2ADA 7 + A 7 2 

A±= i-A^ ± y (i_ A £>2)2 ■ (19) 

Nonetheless, it is possible to see that the subleading prefactor entering in expression 
l|18fl presents an extra cut for A < AI, where, in the infinite time limit, 

A* = - A±£l . (20) 

Note that this cut exists only for A7 > 0, and implies that the right tail of the large 
deviation function of Q± will present an exponential decay, with a slope different than 
the one predicted by the Legendre Transform of n$. The quantitative details on how 
this extra cut will affect 7r(g) will be given later. Let us now give some comments on 
this issue. The energy balance equation for this particular system can be written as: 

AE(t) = E(t) - E(0) = Q x + Q 2 , (21) 

where Q 2 denotes the work performed by the thermostat at temperature T 2 on the 
particle, and is defined by Eq ©. Physical intuition would suggest that for very long 
times, since the Qi's are extensive in time while the "boundary term" AE is not, the 
latter does not come into play when determining the large deviation of, say, Q\, and 
in particular that Q\ and —Q2 have the same large deviation functions. However, 
in several examples ^] El Q H61 H] it has been noted that, if the boundary term 
distribution has exponentially (or slower) decreasing tails then, even in the infinite 
time limit, it can have a fluctuation of order t, and hence is no more negligible. 
Clearly in our case, since the velocity pdf is Gaussian, the energy has exponential 
tails. If the initial condition vo is fixed, then AE has only a right exponential tail, 
which may affect the right tail of the large deviation function of Q±. This is precisely 
the reason for the presence of our extra cut. Besides, if the initial condition is sampled 
on the stationary velocity pdf, the boundary term AE will have two exponential tails, 
and therefore two extra cuts may appear in the A complex plane. When the initial 
condition is sampled in the stationary state, one finds: 

1 

P(A) = J dv e~^P(AM = e* X (cosh(r]t) + - A(A 7 + AD(X - 1) + ' .(22) 

In the infinite time limit this last expression still presents the AI extra cut for A7 > 0. 
Moreover, a new cut appears for AD > A7/2 and for A > A!J_ with 

1 + 2AD - A7 
1 + AD 

These cuts are schematized in a diagram in figure |2 

The Laplace Transform can be safely inverted using the Legendre Transform only 
in the region in which: 

max(A_, X*_) < A < min(A+, A* ) . (24) 

The Legendre Transform gives hence in this region: 



K = \ , An ' ■ (23) 



1 - AD 2 + 2q(Aj - AD) + J (1 + A7 2 - 2A 7 AL>)(1 + 4<f - AD 2 ) 
nc{q) = 2(1 — AD 2 ) 
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Figure 2. (color online) Diagram for the extra cuts in the real A axis. The 
solid lines separate the regions of the parameter space (A7, AD) for which P(X) 
has different domains of definition in the A-real axis. The dashed line is the line 
where T 1 =T 2 - 



If the bounds of P(A) are (on the real axis) A+ and A_ the above expression is valid 
for any q, and the Fluctuation Relation holds. Besides, if the left (right) cut begins 
at A I (A*j_), the Legendre Transform of /«q(A) is valid only for q < q + = fj,' (\*_) 
(q > q- = A*o(^+))- I n this case the expression of ir(q) outside the interval (g_,g+) 
is a straight line ir(q) = a + 0q, and the coefficients a and (3 can be obtained using 
that both 7r(g) and 7r'(g) must be continuous functions of q (this is not evident for the 
derivative 7r'(g), but it is effectively the case in ^1], where essentially the same model 
has been worked out). For example in the region A7 > and AD > A7/2 one has: 



(26) 



with 





' 717 (?) 


q<q- 


n(q) = < 


7T c (g) 


q- <q<q+ 






q> q+ , 


AD 


1 


q - = ~2 


2(2AL> - A 7 ) ' 9+ 



1 A ° (27) 



2A 7 2 ' 

(l + 2AD-A 1 )(l + AD-2g) 
ni {q) ~ 2(1 + A7J) ' (28) 

(l + A 7 )(l + AP + 2 g ) 
nAq) ~ 2(1 + AD) ' (29) 

The behavior of Tr(q) in the others regions can be obtained in a similar fashion. If an 
extra cut starts at X*_, then for q > q+ the ldf ir(q) is equal to ir r (q). Analogously, if 
an extra cut begins at A^, then for q < g_ one has ir(q) — ni(q). In all the other cases 
Tv(q) = TT c (q). The Fluctuation Relation JSJ is only satisfied in the case for which 
n(q) = TT c (q). The Extended Fluctuation Relation is satisfied in the cases where 




Figure 4. (color online) Three dimensional plot of q* , as a function of the 
system's parameters. In the region in which q* = oo the FR is verified. The solid 
line delimitates the region in which q* = 0, where the EFR fails. 



<7_ < < q+ (see the diagram in Fig[3J). In this case q* is the minimum between 
|<7_| and \q+\, and is plotted in Fig0] Surprisingly it also happens that, in a in a 
given region of the parameters' space (i.e. when A7 < 2 AD — 1/AD) one has that 
< q- < q+. In this case q* — and both the FR and the EFR break down. In 
particular it is found, for small values of q, that: 

n(q) - %(-q) = ni(q) - n(-q) = (q , \q\ < q~ (30) 
with £ = 4/Tq — 2/Ti (an illustration is given in Fig|SJ. 
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Figure 5. Plot of the large deviation function n(q) for AD = 0.75 and 
A7 = —0.5. The (red) dotted line shows 7r c (<j), for which the Fluctuation Relation 
holds. The inset shows the plot of n(q) — tt(— q). The (red) dotted line is a straight 
line of slope e, while the (green) dashed line has a slope £ (see text for details). 

We now turn on some remarks concerning the importance of initial conditions. If 
the two temperatures T\ and T2 are equal, it is clear that the average flux (Qi) = 0, 
with w(q) = n(—q). Surprisingly, if the ldf Tr(q\vo) is measured for A7 > (and 
T\ = T 2 ), with a fixed initial condition v , one would see that Tr(q\vo) 7^ Tr(—q\vo) 
(because of the right extra cut), which would lead to the wrong impression that there 
is a nonzero flux from one thermostat to the other, even if they are at the same 
temperature. Of course, if n(q) is measured sampling the initial conditions in the 
stationary state, a new extra cut appears, restoring the symmetry n(q) = n(—q). This 
can be understood with the following remark. The right tail of Tr(q\vo) is somehow 
dominated by the events with an "energetic final condition" (cf. also Jl])- The bias 
enforced by the fixed initial condition (even in the case where vo — (v) = 0), forbids 
the occurrence of "energetic initial conditions" , which would balance (on average) the 
energetic final conditions, and leads then to the appearance of an unphysical flux. 

Our results answer some issues raised in the recent developments of the theory 
of nonequilibrium fluctuations. In particular, we have provided the confirmation that 
if two quantities, both extensive in time (as e.g. Qi(t) and Q2(t)), differ one from 
the other by a "boundary term" (as AE(t)), whose pdf has tails with exponential 
decay, they do not have the same large deviation function. This also leads to a small 
restriction to the validity of the very general result obtained in [HIEj- Furthermore, 
our results clearly show that, even if the Fluctuation Relation (JjJJ is broken, the largest 
eigenvalue no(X) always displays the symmetry property analogous to ijTi)! . While the 
latter result is obviously mathematically robust, this is also the most relevant physics- 
wise, since for small fluctuations around the average and for small nonequilibrium 
drive, it leads to the well-known fluctuation-dissipation theorem. 
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